The crossover from single file to Fickian diffusion 
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The crossover from single-file diffusion, where the mean-square displacement scales as (a;^) ~ t^ , 
to normal Fickian diffusion, where {x^) ~ t, is studied as a function of channel width for colloidal 
particles. By comparing Brownian dynamics to a hybrid molecular dynamics and mesoscopic simu- 
lation technique, we can study the effect of hydrodynamic interactions on the single file mobility and 
on the crossover to Fickian diffusion for wider channel widths. For disc-like particles with a steep 
interparticle repulsion, the single file mobilities for different particle densities are well described by 
the exactly solvable hard-rod model. This holds both for simulations that include hydrodynamics, 
as well as for those that don't. When the single file constraint is lifted, then for particles of diameter 
a and pipe of width L such that (L — 2a) /a = 5c 1 the particles can be described as hopping 
past one-another in an average time thop- For shorter times t <^ thop the particles still exhibit sub- 
diffusive behaviour, but at longer times t 3> thop, normal Fickian diffusion sets in with an effective 
diffusion constant Dhop ~ l/V^fcop. For the Brownian particles, thop ~ S^^ when 5 <C 1, but when 
hydrodynamic interactions are included, we find a stronger dependence than S^^- We attribute this 
difference to short-range lubrication forces that make it more difficult for particles to hop past each 
other in very narrow channels. 

PACS numbers: 



I. INTRODUCTION 

When particles are confined to channels so narrow that 
mutual passage is excluded, the geometric constraints re- 
strict the particles to a single file and a fixed spatial se- 
quence. For short times the mean-square displacement 
may still take the Fickian form (x^) = 2Dot, with a self- 
diffusion coefficient Dq, but at longer times the motion is 
strongly suppressed by collisions with neighbouring par- 
ticles, leading to an asymptotic scaling of the form: 



(1) 



first derived by Harris 1] in a more mathematical con- 
text, and independently by Levitt [2] for particles diffus- 
ing in a narrow pore. Here F is the single file diffusion 
(SFD) mobility. 

There has been an increasing interest in transport 
through highly confined pores [3|, stimulated in part 
by biological realisations such as ion channels Q, and 
aquaporins [^]. Water under extreme nanoscale con- 
finement exhibits behaviour that differs markedly from 
the bulk 0. Single file flow of water is also impor- 
tant in artificial materials like carbon nanotubes Q . 
Transport of simple molecules through porous materi- 
als such as zeolites also show single file sub-diffusive be- 
haviour ^, ilfl, illi] . 

Although the biological and synthetic nanoscale sys- 
tems described above show signatures of SFD, their in- 
terpretation is complicated by numerous other factors 
such as the interaction of the particles with the walls 
of the confining pore. By contrast well defined model 



systems can be created with micron sized colloidal par- 
ticles. One of the major advantages is that the par- 
ticles can be directly imaged in real time with digital 
video microscopy. By using lithography [H, [H, [lj| or 
optical tweezers [l5| to create the one-dimensional con- 
finement for colloidal particles, unambiguous evidence of 
asymptotic SFD {x^) ^ scaling was observed. Lin et 
al. [14] measured the SFD mobility F for different one- 
dimensional packing fractions — pa, where the density 
p — N/Lp, the number of particles is N, the length of 
the pipe is Lp, and a is the colloidal hard-sphere radius. 
They found good agreement between their measured F 
and the SFD mobility for a hard-rod fluid (also known 
as a Tonks gas) (1|: 
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where Ic is defined as the average inter-particle separation 
and tc — l^/2Do is the average time between collisions. 
On time scales t much less than the collision time tc, 
one expects ordinary Fickian diffusion, whereas for time 
scales t ^ tc one expects to observe the asymptotic SFD 
diffusion of eq. H]). 

An obvious question raised by the experiments on SFD 
is what happens as the confinement becomes less severe. 
At some point the system should cross over to ordinary 
Fickian diffusion at long time-scales. This problem was 
first studied by coupling two lattice gas models so that 
particles could jump between chains [16]. Allowing the 
jumps produced a crossover in the mean-square displace- 
ment from subdiffusive \/i scaling to diffusive (cc^) ~ t 
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asymptotic scaling at long times. Difficulties in inter- 
preting measurements of SFD transport in zeolites also 
inspired simulations that exhibited a similar long-time 



crossover to Fickian diffusion [l7j, IS | . 

As shown by Mon and Percus 19| it is convenient to 
analyze this crossover in terms of a Markov process with a 
hopping time thop that measures the average time for two 
particles to pass one another. If t^op S> tc, so that the 
system has reached the SFD regime between hops, then 

the average mean-square displacement of the particles 

1 

scales as If^^p ^ Ftt 



"-hop between hops. Within this picture 
a particle makes on average t/thop hops in a time t so 
that its mean-square displacement scales as 



(^^) ^ ^hop 
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which defines an effective Fickian diffusion coefficient of 
the form 
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The larger thop, the smaller Dhop, since particles passing 
events become more rare. Nevertheless, on time-scales 
t ^ tfiop the final asymptotic scaling will still be Fickian. 
On the other hand, if t^op decreases to the point where it 
is of the order of the collision time tc, then the SFD pic- 
ture is expected to break down and Dhop will approach 
the self-diffusion coefficient Dq. So for the hopping pic- 
ture to be useful, we require t^op ^ tc- 

The hopping time thop has been studied in some depth 
for the case of two discs diffusing between hard walls a 
distance L apart. Mon et al. [l^ [l^, Hf], found a 
scaling of the form 
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as the pore width L approaches the limit where two par- 
ticles can no longer pass. An exponent i> = 2 was found 
both for Molecular dynamics ^2Q] and Brownian dynam- 
ics (B D) l2ll| , but the interpretation of the simulations is 
subtle \22\ . This exponent agrees with direct calculations 
of the diffusion equation l23| , and a simpler transition 
state theory (TST) [13, HjT but not with the effective 
one-dimensional Fick- Jacobs equation [13, [H, [1^ , which 
predicts that — ^■ 

Given that hydrodynamic interactions (HI) can have 
subtle effects on the dynamics of colloidal particles [l^, 
[27| . one might expect there to be differences between 
SFD behaviour measured for colloidal suspensions, and 
the SFD behaviour of particles where hydrodynamics 
is not imp ortant. A careful theoretical study by KoU- 
mann [2a| predicted that the HI do not change the 
asymptotic scaling of Eq. ([1]) , and moreover that the mo- 
bility F can be connected to the short-time collective dif- 
fusion coefficient of the system. Indeed, the experiments 
with colloidal particles cited above exhibit SFD sub- 
diffusive behaviour, as expected. Nevertheless, it would 
be interesting to see if other aspects of SFD behaviour are 



sensitive to HI. Moreover, the subtle crossover from SFD 
to Fickian long-time diffusion for narrow pipes where par- 
ticles can hop past each other, but where thop ^ ^c, could 
well be affected by HI. 



To study these questions, we employ a hybrid MD 
and stochastic rotations dynamics (SRD) computer sim- 
ulation method that has been shown to accurately re- 
produce Brownian and hydrodynamic behaviour for col- 
loidal suspensions Hol- SRD, also known as multi- 
particle collision dynamics, was first described by Male- 
vanets and Kapral [3l| . SRD has been applied to a wide 
number of different systems, including fluid vesicles in 
shear flow [s^l, clay-like colloids [11], sedimentation of 
colloids [34, 35] , colloidal rods in shear flo w I3q| , knots in 
viral DNA [37] and many other examples [38]. We have 
recently used this method to study the role of confine- 
ment on two-dimensional diffusion [s^] and the role of 
finite sized particles on Taylor dispersion (40| . 



To capture the correct asymptotic trends for the 
mean-square displacement, very long simulation runs are 
needed, which explains why so many of the simulations 
in the literature have been on two-dimensional mod- 
els, which are generally faster to simulate than three- 
dimensional systems. In this paper we study strongly 
repulsive colloidal discs in confinement, using both the 
hybrid MD-SRD method as well as simpler BD to com- 
pare what happens when HI are ignored. We find, in 
agreement with experiments on colloidal suspensions [ij] , 
that the simple hard-rod model for the SFD mobility ([2]) 
provides a good fit both for BD and for the simulations 
that include HI. 



When the pipe diameter is such that S = L — 2a > 0, so 
that the mutual passage constraint is lifted, we observe 
a crossover from SFD sub-diffusion to Fickian diffusion 
at longer times t > thop- We measure the distribution of 
hopping times thop, showing that they follow Poissonian 
statistics for simulations with and without hydrodynamic 
interactions. For both models we also find a scaling con- 
sistent with Dhop oc l/y/tiiop- For the BD simulations, 
we find that thop ^ for small S/cr, but thop depends 
more strongly on l/d when hydrodynamic interactions 
are included. We attribute this to repulsive lubrication 
forces, induced when particles approach each other other 
or the walls of the container. 



The paper is organised as follows: In section II, we 
describe the model, our hybrid MD-SRD method and our 
implementation of BD. In section HI we study the SFD 
behaviour of colloids confined to one-dimensional flows 
for different packing fractions 77. Section IV considers the 
case when the mutual passage constraint is lifted, so that 
at long times Fickian diffusion is recovered. Finally, in 
section V, we discuss the main conclusions of this study. 
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II. MODEL AND SIMULATION METHODS 
A. Model 

For the colloid-colloid interactions, we use a very 
strongly repulsive potential of the form: 

^^^^^ " 1 (r > 2V24,) 

where a is the colloid-colloid diameter. The colloid-wall 
interaction is taken to be strongly repulsive as well, and 
a stick-boundary condition is applied to the colloids, as 
well as to the fluid particle-wall interaction in the SRD 
simulations. We use periodic boundary conditions in the 
direction parallel to the pipe walls. 

B. Hybrid MD and SRD simulation method 

To describe the hydrodynamic behaviour of colloids 
induced by a background fluid of much smaller con- 
stituents, some form of coarse-graining is required. The 
hydrodynamics can be described by the Navier Stokes 
equations that coarse-grain the fluid within a continuum 
description. The downside of going directly through this 
route is that every time the colloids move, the boundary 
conditions on the differential equations change, making 
them computationally expensive to solve. 

An alternative to computing a direct solution of the 
Navier Stokes equations is to use particle based tech- 
niques that exploit the fact that only a few conditions, 
such as (local) energy and momentum conservation, need 
to be satisfied to allow the correct (thermo) hydrody- 
namics to emerge in the continuum limit. Simple parti- 
cle collision rules, easily amenable to efficient computer 
simulation, can therefore be used. Boundary conditions 
(such as those imposed by colloids in suspension) are eas- 
ily implemented as external fields. 

In this paper we implement the SRD method first de- 
rived by Malevanets and Kapral [3l|. An SRD fiuid is 
modelled by N point particles of mass m, with positions 
r, and velocities v,. The coarse graining procedure con- 
sists of two steps, streaming and collision. During the 
streaming step, the positions of the fluid particles are 
updated via 



(6) 



In the collision step, the particles are split up into cells 
with sides of length oq, and their velocities are rotated 
around an angle a with respect to the cell centre of mass 
velocity. 



^r,{t + 5ts)= v,.„,,(0 + n,{a) [v,{t) 



.At)] (7) 



where 'Vc.m.i — J^]' (™Vj)/^^-m is the centre of mass 
velocity of the particles in each cell, TZi{a) is the rota- 
tional matrix and St^ is the interval between collisions. 



The purpose of this collision step is to transfer momen- 
tum between the fluid particles while conserving the en- 
ergy and momentum of each cell. 

The fluid particles only interact with one another 
through the collision procedure. Direct interactions 
between the solvent particles are not taken into ac- 
count, so that the algorithm scales as 0{N) with par- 
ticle number. The carefully constructed rotation pro- 
cedure can be viewed as a coarse-graining of particle 
collisions over space and time. Mass, energy and mo- 
mentum are conserved locally, so that on large enough 
length-scales the correct Navier Stokes hydrodynamics 
naturally emerges (sij . It is important to remember that 
for all these particle based methods, the particles should 
not be viewed as composite supramolecular fluid units, 
but rather as coarse-grained Navier Stokes solvers (with 
noise in the case of SRD) [3^. Another advantage of 
SRD is that transport coefficients have been analytically 
calculated 41, 42, 43], greatly facilitating its use. 

The SRD fluid particles can easily be coupled to a so- 
lute as flrst shown by Malevanets and Kapral [1^, and 
studied in detail in a recent paper (30j . To simulate the 
behaviour of colloids of mass M, we use the colloid-colloid 
interaction deflned above, while the solvent particles in- 
teract with the colloids via an interaction of the form: 



h^cs{m''-{^r+i) (r<2vv.) 

1 (r > 2i/Ve.) 



where CTcs is the colloid-solvent collision diameter. We 
propagate the ensuing equations of motion with a Veloc- 
ity Verlet algorithm [46| using a molecular dynamic time 
step At 



R^{t + At) ^ R,{t) + Vi{t)At 



2M 



(8) 



V.(t^At)^V.it)^m±I^^,, (9) 

where Ri and Vi are the position and velocity of the 
colloid, and Fi is the total force exerted on the colloid. 
Coupling the colloids in this way leads to slip boundary 
conditions. Stick boundary conditions can also be im- 
plemented p^ . but for qualitative behaviour, we don't 
expect there to be important differences. In parallel the 
velocities and positions of the SRD particles are streamed 
in the external potential given by the colloids and the ex- 
ternal walls and updated with the SRD rotation-collision 
step every time-step St^. 

The larger the ratio cr/ao, the more accurately the hy- 
drodynamic flow fields will be reproduced. In ref [30j, 
it is shown that using a/aQ — 4.3, and acs — 2ao repro- 
duces the flow flelds with small relative errors for a single 
sphere in a 3d flow. Because we are interested in pro- 
cesses where the colloids can just barely pass each other, 
we chose a finer grid (tr/oo = 8.6). This should in par- 
ticular enhance the resolution of lubrication forces [30j . 
Other parameter choices taken from Refs [30l . [39'| include 
fee = £cs = l.hksT for the colloids, an SRD particle 
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density of 7 = 5 particles per Og and a rotation angle of 
a = \tt. The time-steps for the MD and SRD step are 



set by different physics [30|, and we chose At 
and 5ts — O.lio for SRD, where = ap 
of time in our simulations. 



= 0.025^0 
is the unit 



C. Brownian Dynamics simulation method 

The particle motion for the colloids can also be solved 
using BD. The positions are updated via the equation of 
motion: 

r{t + MBD) = r{t) + ^^F{t) + 5r'', (10) 

where F{t) is the force acting on the colloids, which arises 
from the colloid-colloid interaction, as well as the colloid- 
wall interactions, and each component 6r'^ is chosen from 
a Gaussian distribution with zero mean and variance 
{{5r°f) = 2DoM where Do is the single parti- 

cle diffusion coefficient. With BD, the particles execute 
a random walk with a physical time-scale 



tD 



2Do 



(11) 



called the diffusion time that measures how long it takes 
the particle to diffuse over its diameter. The BD time- 
step was chosen to be AtBD = 2 x lO'^tn, slightly 
smaller than the value used, for example, in a careful 
study by Lodge and Heyes [4^ . We use this conservative 
measure of the time-step because there are indications 
that scaling of the hopping time of colloidal particles 
simulated with BD is sensitive to the time-step [1^; we 
checked that our time-step is small enough. In contrast 
to the SRD method, which also exhibits diffusive motion, 
the BD method does not include HI. 



III. SINGLE FILE DIFFUSION 

We first investigate the SFD behaviour of colloidal 
discs in pipes narrow enough that they cannot pass each 
other. At times less than the collision time tc we expect 
the particles to behave in a diffusive manner, and for 
times t ^ tc we expect SFD sub-diffusion. Lin et al. [T^ 
showed that the following ansatz: 
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2Dot 2Fti/2 



(12) 



provides a good approximation to the crossover from 
Fickian to SFD sub-diffusion. Solving for the mean- 
square displacement gives: 



2Dot 



2Dot 



\ + {Do/F)ty^ l + (i/t^)i/2' 



(13) 




FIG. 1: Log-log plot of the mean square displacement of col- 
loids calculated with BD for a pipe of width L = IAg. Re- 
sults are shown for simulations with increasing line packing 
fractions r] = 0.1, 0.3, 0.4, 0.5, 0.6, 0.7. The straight lines have 
slope 1 and 1/2 respectively and serve as a guide to the eye. 
At long times the colloids show clear signatures of SFD sub- 
diffusive scaling. 
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FIG. 2: Linear plot of the mean square displacement of col- 
loids as a function of ^t/tn, calculated with BD for a pipe 
of width L = 1.4a". Results are shown for simulations with 
increasing line packing fractions 77 = 0.1, 0.3, 0.4, 0.5, 0.6, 0.7. 
The solid lines are a fit from Eq. (|13p . Since Do is fixed by 
the BD simulation parameters the only fit parameter is the 
SFD mobility F. 



by Eq. this crossover time scales as: 



t.-tnl[^)\ (15) 




FIG. 3: Linear plot of tlie m ean square displacement of col- 
loids as a function of ^/TJId, calculated with SRD for a pipe 
of width L = \Aa. In contrast to the BD simulations, these 
include HI. The solid lines are the fit to Eq. p3p . 




FIG. 4: Comparison of the SFD particle mobility extracted 
from BD and SRD simulations, to the for a hard-rod 

fluid from Eq. ([2]). Inset: data at higher packing fractions. 
The results suggest that, within simulation errors, the inclu- 
sion of HI interactions does not strongly afl'ect the ID mobility 
F. 



where we have defined the crossover time 

- (FiDof, (14) 

a measure of the time needed for the system to transition 
from Fickian diffusion to the asymptotic SFD regime. For 
particles that behave like hard rods with a mobility given 



and is connected to the collision time by tx — Itcji^- 
That the two times are essentially the same is not sur- 
prising, given that the (a;^) ~ \Jt behaviour is generated 
by the collisions with neighbours and the long-ranged cor- 
relations these induce. The higher the one-dimensional 
packing fraction 77, the lower t^., and so the more rapidly 
the system should transition to asymptotic SFD sub- 
diffusion. For very small packing fractions the crossover 
time scales with tx/to ~ l/*?^! and so very long simula- 
tions are necessary to observe asymptotic SFD scaling. 
For example, for ry = 0.5, tx ~ 0.63iD, for 77 = 0.1, 
tx ~ 51.5ir, and for r? = 0.01, tx « 6240tr,. 

Simulations were carried out with BD for N = 200 
colloidal particles diffusing between parallel plates a dis- 
tance L — 1.40(7 apart. The one-dimensional packing 
fraction was varied in the range 77 — jf- = 0.1 — 0.7 and, 
after equilibration, simulations were run for total times 
ranging between lOOtu up to 300<d to gather data on 
the mean-square displacement. 

In Fig. [1] we plot the mean square displacement of par- 
ticles undergoing single file diffusion at various packing 
fractions. With the exception of the data for 77 = 0.1, the 
crossover time tx <C lOto so that SFD (x^) ~ \/t scal- 
ing is expected for the entire range of data plotted. The 
measurements are consistent with this scaling. The same 
data is plotted on a linear plot in Fig.[2l where it is com- 
pared to the ansatz of Eq. . Since Dq is given by the 
simulation parameters, there is only one fit parameter, 
namely the SFD mobility F. The plots demonstrate the 
\/t scaling for the mean square displacement v.s. time. 

Simulations were also carried out with SRD. The ex- 
plicit inclusion of a background fluid induces HI between 
the particles, and between the particles and the walls. 
The combination of strong confinement and HI has an 
important qualitative effect on the velocity autocorrela- 
tion functions J39l . and on the Fickian diffusion co- 
efficient 148, 491 of individual particles, especially in two 
dimensions [39l.l5ll|. where HI are particularly important. 
But once the system is confined so strongly that particles 
can no longer pass, the long-ranged HI interactions can 
be strongly screened [5p| . That raises the question: what 
is the effect of including HI on the SFD mobility F? The 
SRD method allows us to address these issues. 

We performed simulation runs for N = 200 colloidal 
particles in a SRD fluid between parallel plates a dis- 
tance L = 1.40(7 and L = 1.98(7 apart. The longitudinal 
packing fractions were the same as for the BD simula- 
tions. After equilibration, simulation data was gathered 
for runs of approximately 10 — SOto in length. Because 
the SRD method includes many solvent particles (0(10^) 
here) , it is computationally more expensive than the BD 
simulations, and so the runs are shorter than those per- 
formed for BD. The mean-square displacements for the 
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FIG. 5: Schematic depiction of the clearance 5 available to 
particles attempting to switch places with one another. Par- 
ticles can hop to pass their neighbours if 5 > 0. 



IV. CROSSOVER FROM SFD TO FICKIAN 
DIFFUSION 

Having worked out some properties of SFD when parti- 
cles cannot pass each other, we now focus on the dynam- 
ics that emerge when the single file constraint is lifted. 
Fig. [5] illustrates how the parameter 5c = S/a = L/a — 2 
describes the maximal distance between the particles 
when they can pass each other. In our case we don't 
quite have hard sphere particles, but the repulsion 
is hard enough that the differences are expected to be 
very small. So we will still use the parameter 5c to de- 
note how easy it is for particles to pass. 



L = lAa case are shown in in Fig. [31 and (a;^) ~ ^/t be- 
haviour is still observed. In this case a fit of the data to 
Eq. (fT3)) was used to extract both Dg and F. The values 
of Dq arc consistent with those of ref. [s^ . 

The values of the fits to F for both BD and SRD are 
shown in Fig. |4] as a function of the packing fraction rj. 
These are compared to the hard rod F^^ of Eq. ^ . At 
all but the lowest packing fractions, good agreement is 
found, suggesting that the inclusion of HI does not signifi- 
cantly affect the value of F/Dq. Moreover, in contrast to 
what is found in two dimensions where HI interactions 
have an important impact on the value of Do [S^, IsH, 
here the measured value of the normalised SFD mobil- 
ity F/ {\/DQa) does not seem to depend significantly on 
the confinement, since differences between what is mea- 
sured at L = 1.4(7 and a.t L — 1.98a appear to be small, 
even though Dq changes considerably [39]. Note that F 
itself does scale with ^/Dq, so the changes in Dq due to 
confinement translate into changes in F, but there are no 
further effects of the HI. Finally, we should point out that 
some further changes in F are expected when L moves 
away from the pure hard rod limit L = a, but these are 
expected to be small [s^ . 

It is only for the lowest packing fraction ?/ — 0.1 that 
we observe significant differences with F^^ for the SRD 
simulations. Interestingly, this deviation is very similar 
to that found by Lin et al. 14] in their experiments. Un- 
fortunately, at this lower packing fraction very long sim- 
ulations are needed to reach the asymptotic SFD regime, 
and so we don't believe that enough data has been gath- 
ered here to reliably extract F. The error bars shown in 
the graph are from a standard non-linear fitting proce- 
dure, but it is very likely that the real error is higher. A 
similar conclusion about the errors was made by Lin et 
al. [13] for their data, which also extended out to t fa 
for the lowest packing fraction. To resolve this question, 
the simulations and the experiments should be performed 
for at least an order of magnitude longer time than they 
have been. 



A. Scaling of the mean-square displacement 

As shown in the introduction, it is very useful to con- 
sider the hopping time thop that it takes a particle to 
switch order with one of its neighbours. From Eq. (g]) it 

follows that the effective Fickian diffusion coefficient is 

_i 

expected to scale as Dhop ~ Ft hop- 
To test this concept we performed a number of simula- 
tions at different values of 5c- As an illustrative example, 
we show the mean-square displacement of = 200 Brow- 
nian particles at 77 = 0.5 and 5c ~ 0.093 in Fig. [SI In this 
case the hopping time thop ~ 86iD ^ tx ~ Q.12tB, so 
that the crossover from the initial Fickian diffusion with 
a diffusion coefficient of Dq to SFD a.i t k, t^ is at much 
shorter times than an eventual longer time crossover back 
to diffusive motion due to the hopping. We simulate for 
up to t = 1040io so that there are on average about 11.4 
hopping processes per particle. At times on the order 
of thop we observe a crossover from {x^) ^ ^/t to the 
asymptotic linear Fickian diffusive scaling (x^) ~ 2Dhopt 
with an effective diffusion coefficient Dhop ~ 0.095DO. 
Although in this case a hopping event occurs only about 
once in every thousand collisions, the net effect is to gen- 
erate a linear diffusion that at longer times leads to a 
significantly larger mean square displacement than if the 
particles couldn't pass. 

Fig. [7| shows a similar set of simulations at 77 — 0.5, but 
now with the SRD methods that includes HI. Although 
the runs are shorter (SRD is computationally more ex- 
pensive), they nevertheless show the same trends as the 
BD simulations, with a crossover to Fickian diffusive mo- 
tion at i > thop with an effective diffusion coefficient 
Dhop ~ O.OSTDq. This scaling behaviour is especially 
clear in the log plot on the top right. For this set of pa- 
rameters, the SRD and the BD give similar results for 
thop and Dhop within the expected errors of the simula- 
tions. 



B. Poissonian statistics of hopping events 

The order of the particles was tracked and the number 
of times that a particle switched places with its neighbour 
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FIG. 6: Crossover from SFD to Fickian diffusion. The top two plots show the mean square displacement of the particles 
calculated with a BD simulation at = 0.5 for 5c = 0.093. The mean hopping time was measured to be thop ~ SGto- On 
timescales of the order of the hopping time, we observe a crossover to Fickian diffusion (x^) ~ t with an effective diffusion 
coefficient Dhop ~ 0.095-Do- This crossover is especially clear in the log plot on the top right. . The bottom two plots show 
the probability Phop that a particle has performed A*' hops in the time intervals intervals, t ~ SAbthop and t ~ llASthop 
respectively. 



was recorded. This was done on an interval large enough 
to avoid double counting short switches during the hop- 
ping process. If the hopping events are independent, then 
the probability that a particle makes N hops is propor- 
tional to the time interval over which a measurement is 
made and should follow Poissonian statistics. The prob- 
ability that such an event will occur n times in a given 
time interval can be expressed as 



P(n;A) = ^^. (16) 



where A = t/thop in our case. The bottom two plots 
of Figs. [6] and [7| show the distributions at two different 
values of A. The solid lines denote the distributions cal- 
culated from the Poisson distribution , and the circles 
denote those obtained from the simulation. Within the 
simulation error, these show good agreement for both the 
BD and the SRD simulations, suggesting that the hop- 
ping events are indeed independent and not correlated 
with each other. 



C. Scaling of thop and Dhop with plate separation 

Further simulations were performed for 1000 Brownian 
particles diffusing in pipes of various sizes. In Fig. [8] we 
plot the hopping times thop and effective diffusion coeffi- 
cients Dhop for BD simulations at 77 = 0.7. The top panel 
shows how hopping time varies with the distance 5c and is 
consistent with thop (5^7^. This exponent agrees with di- 
rect calculations of the diffusion equation [23], and a sim- 
pler transition state theory [13, Hj] , but not with the ef- 
fective one-dimensional Fick- Jacobs equation (20I, [23. [25j . 
Nevertheless, we don't consider our simulations to be on 
a large enough range of Sc to conclusively determine the 
scaling behaviour. Such simulations become increasingly 
difficult for smaller 6c because the hopping events become 
more rare. However, recent Monte Carlo simulations by 
Mon [22] for just two discs are consistent with transition 
state theory over a sufficiently larger range of Sc to con- 
firm the scaling law for that special geometry. His results 
suggest that if full simulations of many discs were per- 
formed over a wider range of Sc they would confirm the 
scaling law we observe over a limited range of Sc- 

In the bottom panel in Fig. [SJ we plot the effective 
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FIG. 7: Crossover from SFD to Fickian diffusion for simulations that include hydrodynamic interactions. The top two plots 
show the mean square displacement of the particles calculated with an SRD simulation ai rj — 0.5 for Sc — 0.093. The mean 
hopping time was measured to be t^op ~ TSto- The bottom two plots show the probability Phop that a particle has performed 
A'' hops in the time intervals intervals, t « 0.8thop(6.3t£)) and t ~ 2.2thop{i7tD) respectively. The straight line shows the 
predicted Poissonian distribution (|16p and the circles are the results from the simulation. 



diffusion coefficient Dhop versus thop- The solid line has 
the slope —1/2 which is consistent with the scaling of 
Eq. (ID), as postulated by Mon and Percus [l^. Again, the 
range of 6c s is not large enough to firmly fix the scaling 
behaviour, but it is at least consistent with D oc t~2 
scaling, and inconsistent with an earher prediction of 
scaling [13] ■ As the pipes become wider and thop becomes 
smaller, we expect Dhop to eventually increase to Dq, the 
self-diffusion coefficient. Thus at larger Sc this simple 
scaling law should no longer hold. 

Simulations were also performed using the SRD 
method to test the effect, if any, of hydrodynamics on 
the hopping process. We simulated a total of 1000 freely 
diffusing colloids in 2d pipes with widths characterized 
by Sc ~ 0.035 — 0.093 The longitudinal colloid density 
was rj — 0.5 in all cases. A separate set of BD simula- 
tions was also performed so that the effect of HI could 
be clearly contrasted and compared. 

The top panel of Fig. [9] shows that the BD exhibits the 
same scaling that was found in Fig. [5] for a higher 
packing fraction. For the larger values of Sc shown, the 
BD and SRD simulations give almost the same value of 
thop/to- However, for smaller Sc, the SRD simulations 
find a significantly higher value of thop, which suggests 



that the HI suppress hopping events. We attribute this 
to the following processes. When two particles come close 
together, or close to a wall, the solvent must be displaced 
for the particle to move past. At very short distances, 
this gives rise to so-called lubrication forces [5^. These 
are repulsive for two particles approaching each other, 
and attractive when two particles are close together and 
then move away from each other (because now solvent 
has to ffow into the space between them). SRD simula- 
tions reproduce the lubrication forces, even on distances 
considerably less than gq [s^ . We thus attribute the in- 
crease of thop with respect to the Brownian simulations 
to lubrication forces, since they make it harder for two 
particles that diffuse towards each other to pass when Sc 
is very small. 

The bottom panel of Fig. [5] shows the effective diffu- 
sion coefficient plotted against the hopping time for the 
four largest separations. Within simulation errors, we 

— 1/2 

recover the Dhop ~ tf^^^ scaling that was also observed 
for the BD simulations. Although the scaling law will 
break down for larger values of Sc, Figs. [8] and [9] sug- 
gest that Dhop will reach a value close to Dq well before 
Sc — 1. The exact value of Sc where Dhop fully converges 
to that of a bulk system will of course depend on packing 
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fractions and other details of the system. 
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FIG. 8: BD simulations of the hopping time tho^ and the 
effective diffusion coefficient Dho^ for packing fraction r\ — 
0.7 and separations Sc = 0.023, 0.035, 0.047, 0.058, 0.070, and 
0.081. Top: Average hopping time plotted versus 1/(5^ The 
straight line is a fit to the data. Bottom: Effective diffusion 
coefficient D^op v.s. t^op The straight line has a slope —1/2 
and serves as a guide to the eye. 
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FIG. 9: Top graph: SRD and BD simulations of the hop- 
ping time thop for packing fraction i] = 0.5 and separations 
5c = 0.035,0.047,0.058,0.070,0.081, and 0.093. For larger 
Sc both simulation methods give the same value of thop, but 
at smaller separations the simulations with HI have larger 
thop's. Bottom: Effective diffusion coefficient Dhop v.s. thop 
from the SRD simulations. The straight line has a slope —1/2 
and serves as a guide to the eye; The data is s consistent with 
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V. CONCLUSION 

We have carried out computer simulations to investi- 
gate single file diffusion and the crossover to Fickian dif- 
fusion when the particle passing constraint is lifted. By 
comparing BD and SRD simulations, we can study the 
effect of HI on these processes. Whereas HI have an im- 
portant effect on the single particle diffusion coefficient 
under confinement, especially in 2d [soj . given Dq, we 
were unable to measure any further effects on the SFD 
mobility F for particles that cannot pass each other. The 
value we measure for F at different packing fractions ttJs 
consistent with the exactly solvable hard rod model 
Experiments on colloidal suspensions [l^, [H, [13] j ^-Iso 
find {x^) fa 2Ft^ scaling and in one case quanti- 
tative agreement with F from the hard-rod model, so 
simulations and experiment both suggest that HI have 
at best a very small effect on the SFD mobility. 

When the no-passing constraint is lifted, we can mea- 
sure the average hopping time thop that it takes a par- 
ticle to switch order with a neighbour. While at shorter 
times tx t <^ thop the particles still exhibit SFD, 
on time scales t > thop, Fickian diffusion emerges with 
{x^) K, 2Dhopt- The effective hopping diffusion coefficient 

scales as Dhop 

We find that for the Brownian simulations thop ^ 
for small (5c, as predicted by TST, but when HI are in- 
cluded, a stronger dependence becomes evident at small 
8c- We attribute this increase of thop to the influence of 
hydrodynamic lubrication forces that make it harder for 
the particles to pass one another. 

The hopping time thop and diffusion coefficient Dhop 
depend strongly on be- This means that measuring how 
thop or Dhop depend on 8^ may be difficult for experi- 
ments, because very small relative errors in the channel 
width can lead to large effects. But this sensitivity to 
changes in be could also be used as an advantage, since it 
could be exploited, for example, by microfluidic applica- 
tions [5^ as a very sensitive way to measure of the size 
of particles in artificial channels. 
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